Photon propagation in a discrete fiber network: An interplay of coherence and losses 
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We study light propagation in a photonic system that shows stepwise evolution in a discretized 
environment. It resembles a discrete-time version of photonic waveguide arrays or quantum walks. 
By introducing controlled photon losses to our experimental setup, we observe unexpected effects 
like sub-exponential energy decay and formation of complex fractal patterns. This demonstrates 
that the interplay of linear losses, discreteness and energy gradients leads to genuinely new coherent 
phenomena in classical and quantum optical experiments. Moreover, the influence of decoherence 
is investigated. 



Although present in almost every experiment, the role 
of photon losses in optics is quite elusive. In linear wave 
optics, losses are often ignored, as in many cases they 
can be scaled out. Mostly, losses are considered to be 
a nuisance which do neither offer any benefits nor intro- 
duce genuinely new coherent dynamics. In this Letter, 
we show experimentally that controlled losses generate 
new effects of linear wave evolution in high dimensional 
discrete networks. 

We demonstrate how coherence and discreteness dras- 
tically impact energy decay in such systems, reducing the 
losses to the sub-exponential regime. Furthermore, the 
interplay between coherence, discreteness and losses leads 
to a completely different intensity distribution over the 
output ports of the network. The resulting formation 
of triangular patterns has similarities to the Sierpinski 
sieve. These experimental results demonstrate that the 
discreteness of both evolution and transverse direction 
has an unexpected influence on lossy photon propaga- 
tion. This gives an entirely new flavor to the role of 
linear losses both in classical and quantum optics. 

Our experiment can be seen as a classical realization 
[H-i] of a quantum walk (QW) of single particles 

and resembles a discrete time version of photonic waveg- 
uide arrays (WGA) |T3] . Recently, WGAs have been uti- 
lized to investigate quantum propagation [l^, In 
quantum optics, many important phenomena are gov- 
erned by interference effects, like the Grover search algo- 
rithm 1^ or the QW of single particles. However, com- 
mon WGAs do not allow for a selective management of 
phases and losses during propagation. The transfer of the 
concept of discreteness to the temporal domain as well as 
an additional discretization of the propagation direction 
in our recirculating fiber-loop setup enables this kind of 
manipulation during each propagation step. This allows 
to explore completely new types of coherent dynamics. 

Moreover, the transport of excitations is largely af- 
fected by the environment in the form of decoherence, 
losses and energy gradients. In photosynthesis as a 




FIG. 1. Principal scheme of all-fiber experimental setup (a) 
without deliberately introduced losses (b) to introduce a con- 
trolled loss of photons by measuring all intensity in the |t) 
state after every second 50/50 coupler; SOA: semiconductor 
optical amplifier; PM: phase modulator; PD: photodiode 



prominent example, the initially coherent transport of 
an exciton that is driven by an energy gradient even gets 
accelerated by the noisy environment lB"12i- Although 
surely present, the influence of particle decay has not yet 
been investigated so intensively. It is therefore desirable 
to gain a better understanding of the interplay of losses 
and spatially varying potentials in the dynamics of co- 
herent processes. 
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FIG. 2. Evolution of wave starting in state \D at position n = in the lossless case (Fig. la); (a-d) measurements and 
simulations of ballistic spreading of intensities \u^f in the upper and in the lower loop without phase modulation; 

dynamics are the same as in a discrete time QW (e-h) discrete time Bloch oscillations for linear phase gradient o = ||- 



In the following, we study experimentally how discrete- 
ness together with photon losses and energy gradients af- 
fect the evolution of light governed by wave interference. 

In our experimental setups (Figs. la,b, see [ll] for de- 
tails), the discrete spatial dimension (position space) is 
emulated by time- multiplexing [l^ • In the first case 
(Fig. la), light propagates in two loops of single-mode 
fiber with a length difference 2AL which are connected 
by a 50/50 coupler. Each roundtrip corresponds to a dis- 
crete propagation step m. A coherent light pulse trav- 
eling in the longer (shorter) loop makes a step to the 
right (left) in position space n and at the central cou- 
pler, multi-path interference takes place. The gain of 
each semiconductor optical amplifier in the loops is ad- 
justed to compensate only the signal losses caused by 
light absorption and monitoring. The amplifiers do not 
affect classical wave interference and the system behaves 
as if it was lossless (see [11] )■ However, they allow us 
to realize a considerably larger field spreading than pre- 
viously reported for QWs [19| or in WGAs, with a big 
potential for even a further increase of propagation steps. 
Certainly, their noise would affect the evolution of non- 
classical light states in the system. 

The distribution of the signal between the upper and 
lower loop is represented by two states |t) and ||). At 
step TO = 0, a pulse is fed into the lower loop (state 
ID). The field distribution after m steps is given by the 
amplitudes in the upper and u,™ in the lower loop 
for a respective position n. Fast photodiodes measure 
the intensities and probed by tap couplers. 

Moreover, a phase modulator is inserted into the lower 
loop to induce a phase or energy gradient in position 
space. 



the algebraic equations 
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where a is the induced phase shift between two positions 
n. In addition, a phase shift of ^ acquired in the coupler 
has been taken into account. The discrete time QW of 
single photons is governed by the same equations [I!, 0] . 
Moreover, the dynamics in our system are very similar 
to classical photon evolution in WGAs [H , 2^, 21 1 . 

Without phase modulation (a = 0), we observe a bal- 
listic spreading of the light field (Fig s. 2a-d) as it is well- 
known from continuous WGAs [24I and QWs [4|. The 
initial asymmetry which is caused by injecting the pulse 
only into the lower loop remains conserved (Figs. 2c, d). 

Applying a phase shift (a ^ 0) which grows linearly in 
position to the ||) state at every propagation step to leads 



to discrete-time photonic Bloch oscillations [23[. Similar 
to WGAs [3, U^, the field distribution recovers in a 
quasi-periodic fashion, as it is nicely demonstrated by our 
experiments (Figs. 2 e,g) and confirmed by simulations 
(Figs. 2 f,h). 

So far, on the large scale our lossless system shows 
the same effects as continuous-time WGAs which is also 
expected from the fact that discrete-time and continuous- 
time QWs have the same limiting distribution 2^, 21 1. 

The situation changes completely if we introduce con- 
trolled losses by periodically removing all intensity in the 
It) state. We realize this by a destructive measurement 
of the It) intensity every second time the light has passed 
the central 50/50 coupler. In a slight modification of our 
measurement setup (Fig. lb), we add a second 50/50 



The evolution of the amplitudes can be described by coupler and connect the two couplers with fiber pieces 
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FIG. 3. Evolution of wave starting in state \D at position n = in the case with losses (Fig. lb); (a,b) measured and simulated 
intensities \a^f in the loop with no phase modulation and (c,d) intensities \q. 



simulated fractal patterns of intensity 



at the output port; (e,f) measured and 
for linear phase gradient a — net gain g ^ 3.4/step; (g,h) a — ||-, g k: 3.3 



that have a length difference 2AL, thus realizing the step 
in position space. Only one of the output ports of the 
additional coupler is fed back into the loop, whereas all 
intensity in the second output port is directed to a photo- 
diode. Here, a step m corresponds to a full roundtrip of 
the light in the loop, i.e. to passing both 50/50 couplers 
and the light propagation is governed by the equations 



(2) 



n+lj 



where a™ is the light amplitude inside the loop and g™ 
is the amplitude that is dissipated at the output port at 
position n and step m. Again, we assume that all exper- 
imental losses stemming from absorption or monitoring 
are exactly compensated by the amplifiers. However, a 
compensation of the photons that are intentionally dissi- 
pated at the output port is not possible, because whole 
pulses leave the system and all information which they 
carry (amplitude and phase) is completely lost (see fl8|). 
Due to its transverse degrees of freedom, the loss induced 
evolution in our system is by far more complex than that 
observed in e.g. two coupled harmonic lossy oscillators 

In the absence of phase modulation (a = 0), equation 
(2) is similar to a discretized diffusion equation for the 
amplitudes a™. Moreover, the output amplitudes g™"^^ 
are the discrete derivative of a™. For large iterations 
(m ^ 1), the distribution of amplitudes converges to a 

Gaussian shape a™ « ,'7 exp 



\ 2m )■ 



The ballistic 



spreading in the lossless case is thus transformed into a 
diffusive one. Our measurements (Fig. 3a) demonstrate 
a linear increase of the variance (Fig. 4a) which confirms 
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FIG. 4. Transition to a classical random walk by adding de- 
coherence in the form of random phase shifts; (a) variance of 
Gaussian fit to measured intensity distribution in the coher- 
ent lossy system (Fig. lb) (blue circles) and to the measured 
average distribution in the case of decoherence (red crosses); 
lines: simulation of the coherent system and of a classical 
random walk; (b) fractional energy loss per step to visualize 
the low losses of the coherent system and the loss of 1 /2 per 
propagation step m in the case of decoherence. 



the diffusive evolution. At the output port we observe 
the distribution's discrete derivative (Fig. 3c). Again, 
the agreement with simulations is very good (Figs. 3b, d). 

One remarkable feature of this system is a surprisingly 
low rate of photon losses. Looking at the setup in Fig. 
lb, one could expect that at each roundtrip half of the 
photons are coupled out to the photodiode, correspond- 
ing to a loss of 1/2 per propagation step. But owing to 
wave-mechanical interference, the outcome is completely 
different. Assuming the absence of all losses except for 
the selectively introduced ones at the output port, the 
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energy E{m) J2 I 



that remains inside the loop at 



step m decays as slow as E{m) 



as is found by 



integrating the Gaussian shape. This means that E{m) 
decays slower than any exponential function despite the 
open port of the second 50/50 coupler in Fig. lb. The 
energy Q{m) X) I?™! ^ "T^ ^^at is dissipated at 

the output port at step m decreases much faster than 
E(m), leading to a fractional loss of '^^J?^^'^ ~ — which 
is quantitatively confirmed by our measurements (Fig. 
4b). A system with periodic boundary conditions in posi- 
tion space even possesses a completely lossless mode (see 
(l8|). It is important to note that the photon losses in 
our system do not lead to decoherence [9|, |19|, mHi, but 
instead a high degree of ordering can be observed. The 
astonishing stability of the field in the loop is quite ana- 
log to the quantum mechanical Zeno effect [23|, where a 
permanent measurement prevents a state from decaying. 

If we intentionally remove the coherence from our sys- 
tem by adding random phases [l^ to the light pulses 
at every step and position, a photon gets lost at the out- 
put port with a probability of 1/2 per step as expected 
classically (Fig. 4b). The average dynamics approach a 
classical random walk, where the intensities |a™| per- 
form the diffusion process instead of the amplitudes a™. 
This acceleration due to phase noise 17] is also supported 
by a roughly twofold increase of the variance (Fig. 4a). 

The introduction of a phase gradient a — ^tt, with 
coprime p and q, also leads to an exponential decay of 
E{m). In contrast to the lossless system (Figs. 2e-h), 
no Bloch oscillations are observed. Instead, we discover 
a completely new and unexpected behavior in the case 
of losses: The intensity |a™|^ evolves in a triangular pat- 



tern. The field a™ below the basis of each triangle of 
length q is exactly zero because of symmetries in the in- 
terference of all possible photon paths (see [Si)- With 
increasing p and q, the intensity pattern inside the trian- 
gles acquires a complex fractal nature with similarities to 
the patterns known from linear cellular automata jsB - fs^ . 
An analytical approach to further describe the observed 
patterns will be discussed elsewhere. Our measurements 
shown in Figs. 3e,g and 5a, c are, to our knowledge, the 
first experimental demonstration of this type of fractal 
pattern in optics. The rich structure of the arising frac- 
tals is confirmed by simulations (Figs. 3f,h and 5b, d). 

In summary, we demonstrated an unexpected role of 
photon losses in an optical system that is discrete in 
both time and space, linking our observation to photonic 
waveguide arrays and quantum walks. The introduc- 
tion of state-selective losses leads to the new phenomena 
of sub-exponential losses and fractal evolution patterns 
which have no direct counterpart in continuous WGAs 
or QWs. All discussed effects have been demonstrated 
experimentally and confirmed by simulations. Our work 
gives a new perspective to particle losses in classical and 
quantum systems. This will trigger a thorough search 
for other types of such new effects in linear and nonlin- 
ear optics, possibly with the help of methods known from 
cellular automaton research. Apart from this, the exper- 
imentally realized fiber loop schemes might likewise be 
applied to control multi-pulse regimes in fiber ring lasers. 
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Excellence Engineering of Advanced Materials (EAM). 



5 



[2 

[3; 

[4] 

[6: 

[7] 
[8 

[9 

[lo: 
[11 

[12 
[13 



* |ulf.peschel@mpl.mpg.de| 
P. Knight, E. Roldan, 



and 



Sipe, 



'Phys. Rev. A 68, 020301 (2003) ' 
D. Bouwmeester, I. Marzoli, G. Karman, W. Schleich, 
and J. Woerdman, Phys. Rev. A 61, 013410 (1999) 
P. P. Rohde, A. Schreiber, M. Stefanak, I. Jex, and 

C. Silberhorn, New J. Phys. 13, 013001 (2011) 
J. Kempe, |Contemp. Phy s. 44, 30 7 (2003), 

Y. Aharonov, L. Davidovich, and N. Zagury, 
Fhys. Rev. A 48, 1687 (1993) 

A. Schreiber, K. N. Cassemiro, V. Potocek, A. Gabris, 
P. J. Mosley, E. Andersson, I. Jex, and C. Silberhorn, 
|Phys. Rev. Lett. 104 (2010) , 

M. Karski, L. Forster, J.-M. Choi, A. Steffen, W. A lt, 

D. Meschede, and A. Widera, 'Science 325, 174 (2009) 
H. Schmitz, R. Matjeschk, C. Schneider, J. Glueck- 
ert, M. Enderlein, T. Huber, and T. Schaetz, 
Phys. Rev. Lett. 103 (2009) 

Fedrizzi, B. P. Lanyon, 
and A. G. White, 



"M. X. Broome, A. 
L Kassal, A. Aspuru-Guzik, 



Phys. Rev. Lett. 104 (2010) 



F. Zahringer, G. Kirchmair, R. Gerritsma, E. Solano, 
R. Blatt, and C. F. Roos, Phys. Rev. Lett. 104 (2010) 
D. N. Christodoulides, F. Lederer, and Y. Silberberg, 
Nature 424, 817 (20 03) 

Y. Bromberg, Y. Lahini, R. Morandotti, and Y. Silber- 
berg, 'Phys. R ev. Lett. 102 (2 009) 

A. Peruzzo, M. Lobino, J. C. F. Matthews, N. Matsuda, 
A. Politi, K. Poulios, X.-Q. Zhou, Y. Lahini, N. Ismail, 
K. Worhoff, Y. Bromberg, Y. Silberberg, M. G. Thomp- 
son, and J. L. OBrien, [Science 329, 1500 (2010) | 



[14 

[15; 
[16: 

[17 

[is: 

[19 

[20 

[21 
[22: 

[23: 

[24 

[25: 

[26 
[27 

[28' 
[29 

[30 
[31 
[32: 



P. G. Kwiat, J. R. Mitchell, P. D. D. Schwindt, and 
A. G. White, J. Mod. Opt. 47, 257 (2000) , 
M. B. ' Plenio and F. Huelga, 

'New J. Phys. 10, 113019 (2008)J 

M. Mohseni, P. Rebentrost, S. Lloyd, and A. Aspuru- 
Guzik, J. Chem. Phys. 129, 17411)6 T2OO8) 
F. Caruso, A. W. Chin, A. Datta, S. F. Huelga, and 
M. B. Plenio, J. Chem. P hys. 131, 105106 (2009) 
See supplemental material at [link to be inserted] for a 
detailed explanation of the experimental setup. 
A. Schreiber, K. Cassemiro, V. Potocek, 
A. Gabris, I. Jex, and C. Silberhorn, 

Phys. Rev. Lett. 106, 180403 (2011)j 
F. Strauch, Phys. Rev. A 74 (2006), 



A. M. Childs, Commun. Math. Phy s. 294, 581 (2009) 
T. Pertsch, T. Zentgraf, U. Peschel, A. Brauer, and 
F. Lederer, Phys. Rev. Lett. 88 (2002) | 
A. Wojcik, T. Luczak, P. Kurzynski, A. Grudka, and 
M. Bednarska, Phys. Rev. Lett. 93 (2004) 
R. Morandotti, U. Peschel, J. Aitchison, H. Eisenberg, 
and Y. Silberberg, P hys. Rev. Lett. 83, 4756 (1999) 
T. Pertsch, P. Dannberg, W. Elflein, A. Brauer, and 

F. Lederer, Phys. Rev. Lett. 83, 4752 (1999): 

G. Weinreich, J. Acoust. Soc. Am. 62, 1474 (1977)] 

F. AbduUaev, V. Konotop, and V. Shchesnovich, 
Physical Review A 83 (2011) , 



V. Kendon, Math. Structures Comput. Sci. 17, 1169 (2007) 

M, 



Goniilo l, E. Aydiner 
|Phys. Rev. A 80 (2009)} 



and O. Miistecaplioglu, 



S. Wolfram, Rev. Mod. Phys. 55, 601 (1983) 
S. Takahashi, J. Comput. Syst. Sci. 44, 114 (1992)] 
S. Willson, Discrete Appl. Ma th. 8, 91 (1984]| 



